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'nJ" , Abstract 

The fermion flavor Nj dependence of non-perturbative solutions in the 
strong coupling phase of the gauge theory is reexamined based on the in- 
Qh! terrelation between the inversion method and the Schwinger-Dyson equation 

^ ' approach. Especially we point out that the apparent discrepancy on the value 

of the critical coupling in QED will be resolved by taking into account the 
higher order corrections which inevitably lead to the flavor-dependence. In 
r> I the quenched QED, we conclude that the gauge-independent critical point 

j^ ■ etc = 27r/3 obtained by the inversion method to the lowest order will be re- 

duced to the result Oc = vr/3 of the Schwinger-Dyson equation in the infinite 
order limit, but its convergence is quite slow. This is shown by adding the 
chiral-invariant four-fermion interaction. 
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1 Introduction 

To study the dynamical symmetry breaking, we need some non-perturbative tool. In 
recent several years, the Schwinger-Dyson (SD) equation has played such a role and 
much effort is devoted to solving the SD equation for the fermion propagator in gauge 
theories. In QED4 with Nf fermion flavors, QED4[A^/], there exists a critical gauge 
coupling Oc for a, above which [a > ac) the chiral symmetry is spontaneously broken, 
i.e., {ip'4') 7^ in the absence of bare fermion mass. Such non-perturbative solutions of 
the SD equation ware able to reveal many interesting features of the strong coupling 
phase in gauge theories, see e.g., [|l|, for review. However this approach has the 
following drawbacks. 

1. systematic improvement: It is rather difficult to incorporate the higher order 
effect systematically. In contrast to the perturbation theory, it is not neces- 
sarily clear which terms should be taken account of to the next order. Even 
if this may be possible, it is quite difficult to actually solve it without drastic 
approximations. In essence, the SD equation for the fermion propagator should 
be solved simultaneously with those for the photon propagator and the vertex 
function and such an attempt has been done in a restricted case with a bare 
vertex [Q. However the ad hoc approximation for the vertex function and the 
photon propagator inevitably truncates this series of hierarchy. 

2. gauge invariance: Any gauge-independent result has not been obtained so far. 
Almost all the analyses of the SD equation have been done under the choice 
of the Landau gauge. However, qualitatively different and pathological results 
are obtained, once the other gauges than the Landau gauge are chosen |Q]. 
This implies that the adopted approximation breaks the gauge invariance which 
should hold in gauge theories. In other gauges than the Landau gauge, we must 
take into account the wavefunction renormalization for the fermion, even in the 
quenched approximation ^. This enforces us to consider the vertex correction 
in the light of the Ward-Takahashi (WT) identity. Indeed, many physicists have 
tried to recover the gauge-invariance by modifying the vertex so as to satisfy 
the WT identity. However the WT identity does not uniquely determine the 
vertex function of the SD equation in the non-perturbative sense @ , in contrast 
with the perturbation theory. We must impose further conditions or require- 
ments; absence of the kinematical singularity, multiplicative renormalizability, 
agreement with the perturbation theory in the weak-coupling limit and so on. In 
spite of these trials, however, completely gauge-independent results have not yet 
been obtained in the framework of the SD equation. See references [^ |, |, [10 
for the most recent result. 



The inversion method [|T2| is a generahzation of the Legendre transformation in 
the effective-action formahsm [^ and enables us to calculate the order parameter in 
the symmetry-breaking non-perturbative phase in a perturbative way. The order pa- 
rameter in the effective action is uniquely determined by the form of the source term 
in the action Sj from which we start, while in the inversion method it can be chosen 
more freely providing that it coincides with the original action when the source term 
vanishes. In the inversion method, we can extract the non-perturbative solution by 
perturbative calculations. This method allows us to perform the systematic improve- 
ment [0. By using this method, the strong-coupling phase of QED is studied in a 



gauge-invariant manner |1^, |I^. Accordingly the inversion method is very useful and 
has remarkable features. 

However it has not attracted so much attention of most physicists who have been 
studying the dynamical symmetry breaking. For one thing, the relationship between 
the inversion method and the SD equation approach was not necessarily made clear so 
far; for another thing, the SD equation approach is more direct to obtain the physical 
quantities, such as the fermion mass and the pion decay constant. 

Moreover there exists an apparently contradicting result on the critical coupling 
already in QED between the inversion method and the SD equation approach. 

In this letter we show that there is no contradiction between the result of the 
inversion method and the SD equation approach. Our systematic treatment of the 
solution of the SD equation and its comparison with the inversion result suggests that 
the correct gauge-independent critical coupling in the quenched limit should be given 
by ttc = 7r/3. 

2 SD equation and inversion method 

First we recall the result of the SD equation. The SD equation is also derived from 
the inversion method as follows. After introducing the source term. 



Sj = J d^xj d^yJix,y)^ix)^{y), (1) 

into the original action of (euclidean) QED4[A^j]: 

SoED = d xCc 



QED — / " ■^J-'QED, 

Cqed = -^F,,F^'' + r{^p-mo)r + eori,rA^, {a = l,...,Nf), (2) 



the inversion process with respect to the general bilocal source J{x, y) yields, to the 
lowest order O^Cq), the equation: 

J{p) = -Sipy' + Soip)-' - elj^^i,S{q)^^D>^^{q-p\ (3) 

where J{p) is the Fourier transformation of J{x — y), and Sq{p) and DQ^{k) denote the 
bare fermion propagator and the bare photon propagator, respectively. By putting 
J = 0, the quenched ladder SD equation for the fermion propagator S{j>) is obtained. 



The quenched ladder SD equation in the Landau gauge fl^] leads to the critical 
coupling ^: 

a,(0) = ^ = l. 041975 ■■•. (4) 

On the other hand, the lowest order inversion method applied to the chiral con- 
densates ( TpTp) of QED4[A^j = 1] leads to the gauge- independent critical point 0: 

ac= — = 2.094395..., (5) 

as shown in [|^. The critical point Uc depends in general on the number of fermion 
flavors Nf, adNf) ||l^ [T^, ^. To the lowest order inversion, however, the critical 



coupling ttc of QED4[A^j] is independent of the flavor number of fermion, Nf, as 
pointed out in ||15|| . 

The critical coupling eq. (^ is twice as large as the result eq. (H) of the quenched 
ladder SD equation eq. (^) in the Landau gauge. Then the result eq. (^ seems 
to contradict with the inversion scheme. Because the result eq. (^) follows from a 
special case, i.e., the uniform source J with J{x, y) = J6^{x — y) in the same inversion 
procedure in which case we have the relation after inversion: 



An 



NfAj 



1-^ 

etc. 



(^^) 



^The value adA) of the critical coupling obtained from the quenchd ladder SD equation for the 
fermion propagator S{p) = [A(p^)7^p^ — B{p'^)]~^ depends on both the gauge parameter ^ and the 
ultraviolet cutoff A [Q ^. For ^ < —3 there is no nontrivial solution: B{x) = 0. For > ^ > —3, 
ac(A) exists, but ac{A) | as A t cxo. In the Landau gauge ^ = 0, ac(A) — > 7r/3 as A f oo. For 
^ > 0, ac{A) -^ a^ > 7r/3 as A | cx) and a'^ increases monotonically in ^. Note that A{p^) > 1 
for £, > and < A{p'^) < 1 for i^ < 0. The quenched ladder SD equation with the bare vertex is 
consistent with the WT identity only in the Landau gauge where A{p^) = 1. 

^This value is obtained if we put the same ultraviolet cutoff for fermion momentum Ay and 
photon one Ap in ac = 27r/(3?7) with r] = A^/A^ |14|. Regularization independence of the critical 
coupling will be discussed in a forthcomming paper in detail. 






TiNjA}'^"^' \N]AjAl 
+0{{ilj^)Hn{tlj^f). (6) 

Rather the inversion result eq. (|^) is in good agreement with that of the un- 
quenched ladder SD equation (in the Landau gauge) for one fermion flavor Nf=l 
[0,[T|, 0,0,0: 

ac{l) = 2.00, (7) 

in the presence of the 1-loop vacuum polarization for the photon propagator 

where A is the ultraviolet cutoff. 

How can we make a compromise among these results? Or, different formalism may 
give the different critical point? How can the effect of photon vacuum polarization be 
included and the unquenched result be extrapolated to the quenched hmit, U^^^k'^) -^ 
0, i.e., D^^'ik) -^ D^^ik)? 

We clarify these points based on the systematic treatment of the SD equation 



13,123, |2^ 



3 Analysis 

All the above statements will be made clear if we include the chiral-invariant four- 



fermion interaction [2^| (gauged Nambu-Jona-Lasinio model |^) as 



c,F = ^Gonrrf - iri^rf], {a = i,...,Nf), (9) 

and consider the critical line in the two-dimensional space of coupling constants 
(co. Go) instead of the critical point. We define the dimensionless four-fermion cou- 
pling g hj g = -^^ — , so that the critical four-fermion coupling constant g^ is kept 
fixed: gc = ^ for arbitrary Nf (in the chain approximation). 

In the previous paper ||15[ the inversion method was applied to the gauged Nambu- 
Jona-Lasinio (NJL) model which reduces to QED by switching off the chiral-invariant 
four-fermion interaction. Q To the lowest order, we have obtained the gauge-independent 
results: 



•^Within the framework of the SD equation, an approach to recover the gauge invariance approx- 
imately was tried for the gauged NJL model in Efl] . 
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1. the critical line separating the spontaneous chiral- symmetry breaking phase 
from the chiral symmetric one, 

2. the mean-field value 1/2 for the critical exponent of the chiral order parameter 

3. the large anomalous dimension 7^ = 2 for the composite operator iptp on the 
whole critical line. 

These results are in good agreement with the previous results of the SD equation 
which includes the 1-loop vacuum polarization to the photon propagator in QED 



2^, ^. However the critical line obtained from the SD equation deviates from the 
inversion result in the region of the strong gauge coupling, a = Cq/Att > 4 where the 
bare coupling of the four-fermion interaction is negative [|^. This indicates necessity 



of incorporating the higher orders in the inversion calculation. 

From the SD equation, the critical line of the gauged NJL model is shown to take 
the following form P^l E^]. 



where 






Oi = — O"^ — 1, 



a 



n(>3) 



1 \n n—1 



^' i=2 



t 1 \n n 

bni>2) = ^—^i^-nfUP^'i^)^ (11) 



^! i=2 



and 



with 



P,(a) := a' - 2{i - l)a + {i - l){i - 2), (12) 



9 ^TT 



This is derived from the asymptotic solution of the SD equation [0 with an expansion 
parameter Zq^. Up to the order 0{zq^^), the critical line is written as 



2\ ,9 / 3 2 
-A2-21-- + - 
cr/ V cr cr" 



g = l-2X-{l-^]X^-2{l-- + — )X^ 



5_?H + 52^14^,._2A_2V^-11 + 4)A= 



/ 2\V 187 379 318 90\ ,. 

-4 1-- 33 + + A^ 

V a/ \ a a^ o-^ a^ / 

2\V 3304 9960 14508 9996 2520\ ,o 
1__ 429 + — — + — -)^' 

( 2\2/ 7048 28760 61756 72708 
-2 1-- 715 + ^ ^ + ^^ 

43608 10080\ 9 



a^ a"" 



29479 153724 445060 766920 



2 1-- 2431 + , :— + 



aJ \ a cr^ a^ o"'" 
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776940 420120 90720 \ 
0-5 0-6 a^ y'' 

+0 (\'') , (14) 



where 



A:=- = 5^. (15) 

Now this expression allows us to take the quenched limit Nf —>■ 0: 

g = 1-2A-A2-2A^-5A^-14A^ 

-42A^ - 132A^ - 429A^ - 1430A^ - 4862A^° - C(A"), (16) 

with A being kept fixed, although both a and zq have singular A^j-dependence. There- 
fore we are lead to conclude that the critical line eq. (0) has the correct quenched 
limit g 11: 

^Similar argument can be done for the QCD-like gauged NJL model, as already pointed out in 
l28|,i. 



since eq. (O) is expanded as 



,.1_2A-V!1S(^^A». 



oo 



n=2 



To the lowest order 0{a), the critical line eq. (|l^) is given by 

9=1-^. (19) 

This is independent of the fermion flavor and completely agrees with the lowest order 
inversion result |T^. This implies the critical coupling ac = 27r/3 in pure QED {g = 0) 



which is nothing but the result of Ukita, Komachiya and Fukuda[|14|. 

The complete agreement of the SD result eq. (|19]) with the inversion result in the 
lowest order 0{a) is somewhat surprising. Because in the analytical treatment of 
the SD equation the standard approximation Il^^^{(p — g)^) = n^^)(max{p^, g^}) is 
adopted in order to perform the angular integration and to conclude no wavefunction 
renormalization in the Landau gauge as in the quenched case, while there is no ap- 
proximation in the inversion method to this order. However the critical coupling is 



insensitive to this approximation |]T9l, although numerical calculation of the SD equa- 
tion leads to slightly different value under different approximations [^. For Nf = 1, 
our result shows the critical point up to the order 0{a^^) (see Table 1) 

a^(l) = - X 1.90942 = 1.9995. (20) 

In order to obtain A^j-dependent result, we must take into account at least the 
next order (9(q;^). Actually, the critical line to the order O(a^) is given by 

From the viewpoint of the inversion, the different A^/-dependence does appear in the 
next order from a vacuum diagram with two fermion loops, see Figure 1 of ||14 |. 



It is worth remarking that the numerical result of the SD equation, ac = 2.00 for 
Nf = 1 |]l^ , is obtained from the above equation by solving the second order algebraic 



equation: = 1 — 2A — (1 — |A^/)A , which has indeed a positive solution at ac = 
2.03924 for Nf = 1. The ratio of the second order term to the first order one is 
(1 - |iV/)/4 = 1/36 = 0.0277 for A^^ = 1 at a = 27r/3. This is why the lowest order 
result eq. (^ is very close to the value eq. (|^). Therefore this coincidence for Nf = 1 
is an accidental one due to the negligible second order term. However, the coefficient 
of the C(a^) term in eq. (|2T| ) changes the sign for Nf > 9/8. Hence g y -|-cx) as 



a —>■ +00, which contradicts with the SD equation resuh: g \ —00 as a ^ +00 



[^, |30|. Nevertheless the critical line eq. (|2l|) can give good estimate at least for not 
so large a: a < 0{ac)- Actually the result in the case of Nf = 2 should be compared 
with the value ac(2) = 2.5 obtained from the numerical solution of the SD equation 



n,ig,o 



We must mention the behavior of the series in more higher orders. For Nj > 0, 
it is known ||22|, § that the series eq. (|l^) is not convergent and at most strong 
asymptotic series (hence Borel summable) [^, while in the quenched case the series 
eq. (|16D converges to eq. (p!7|) as N ^ 00. Therefore the critical coupling constant 
ac approaches to the optimal value at a certain N and then deviates from it, which 
can be seen in Table 1. However, for large Nf the asymptotic series fails to give good 
estimate on the critical coupling, since the expansion parameter Zq^ should be small: 



z. 



-1 







Nja/^Sn) < 1, as already pointed out |T^ 



In Table 1, it is demonstrated that all order results are needed in order to recover 
the correct critical coupling adO) = 7r/3 of the quenched limit from the unquenched 
SD equation, since the convergence to ac(0) = vr/S is very slow; adO) = (vr/S) x 
1.00957 to the order 0{a^°°). 

4 Conclusion 

For QED4[A^j], it has been shown that the flavor A^j-independent critical value ac = 
27r/3 is a common value to the "lowest order" f\ 0{a) in the SD equation approach as 
well as the inversion method. The key point lies in the fact that the A^/-dependence 
appears in the next order 0(q;^) for the first time. Therefore it is merely an accidental 
coincidence that the lowest order inversion result ac = 27r/3 is nearly equal to the 
value ac = 2.00 obtained from the unquenched ladder SD equation for Nf = 1. 

The quenched case, Nf = 0, can be extrapolated from the expression eq. (|1^ 
which we have obtained from the SD equation. However, the convergence to adO) = 
7r/3 is very slow as demonstrated in Table 1. Therefore, in order to obtain the 
quenched result eq. (^ in the inversion method, we need to include all the orders of 
vacuum diagrams with only one fermion loop (with no further internal fermion loop). 
This subtlety seems to reflect the essential singularity in the Miransky scaling [jl6| in 
the quenched QED. To obtain such a correct gauge-independent quenched result, we 
have to sum up all the order results of the inversion. 

To really confirm our claims, anyway it is necessary to perform the calculation 
to higher orders in the inversion. It is easily shown that the gauge-independence 
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The " order" of the resuh obtained from the SD equation is defined in correspondence with that 



of the inversion method, see eq. (14) 
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is guaranteed also in the higher orders. Then it is important to study the next 
order O(a^) of the inversion method in order to determine whether or not the gauge- 
independent inversion result really agrees with the SD equation result in the Landau 
gauge eq. (0). Details on this point will be reported in a subsequent paper. 
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order A^ Nf = N 



f 



'f 



Nf 



1 


2.0 


2.0 


2.0 


2 


1.65685 


1.94733 


2.71849 


3 


1.50434 


1.92136 


3.08042 


4 


1.41548 


1.91234 


2.46357 


5 


1.35635 


1.91018 


2.87205 


6 


1.31374 


1.90964 


2.52280 


7 


1.28136 


1.90947 


— 


8 


1.25579 


1.90943 


2.33712 


9 


1.23501 


1.90942 


— 


10 


1.21775 


1.90942 


2.00813 


15 


1.1616 


1.90938 


— 


20 


1.13025 


1.91519 


0.92838 


25 


1.10993 


1.69725 


— 


30 


1.09556 


— 


0.56558 


35 


1.08481 


1.12770 


— 


40 


1.07643 


* 


* 


50 


1.06413 


* 


* 


60 


1.05549 


* 


* 


70 


1.04906 


* 


* 


80 


1.04406 


* 


* 


90 


1.04006 


* 


* 


100 


1.03678 


* 


* 


200 


1.02077 


* 


* 


300 


1.01478 


* 


* 


400 


1.01158 


* 


* 


500 


1.00957 


* 


* 



Table 1. Critical coupling constant Oc in unit of 7r/3 in QED4[A^j], obtained by solving 
the algebraic equation eq. (|T^ g = f{a) = up to the order 0{a'^). To the lowest 
order 0{a), a^ = 27r/3 irrespective of the fermion flavor Nf. In the quenched case 
{Nf = 0) the critical coupling approaches the value, ac = vr/3 for large A^. The 
symbol — stands for no real positive solution to the algebraic equation f{a) = and 
* implies that we have no available data. 
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